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We analyse numerically the linear stability of a liquid metal flow in a rectangular duct 
with perfectly electrically conducting walls subject to a uniform transverse magnetic 
field. A non-standard three dimensional vector stream function/ vorticity formulation 
is used with Chebyshev collocation method to solve the eigenvalue problem for small- 
amplitude perturbations. A relatively weak magnetic field is found to render the flow 
linearly unstable as two weak jets appear close to the centre of the duct at the Hartmann 
number Ha ~ 9.6. In a sufficiently strong magnetic field, the instability following the 
jets becomes confined in the layers of characteristic thickness S ^ HaT^^'^ located at 
the walls parallel to the magnetic field. In this case the instability is determined by 5, 
which results in both the critical Reynolds and wavenumbers numbers scaling as ^ . 
Instability modes can have one of the four different symmetry combinations along and 
(— I ■ across the magnetic field. The most unstable is a pair of modes with an even distribution 

^ ^ of vorticity along the magnetic field. These two modes represent strongly non-uniform 

vortices aligned with the magnetic field, which rotate either in the same or opposite 
senses across the magnetic field. The former enhance while the latter weaken one another 
provided that the magnetic field is not too strong or the walls parallel to the field are 
not too far apart. In a strong magnetic field, when the vortices at the opposite walls are 
well separated by the core flow, the critical Reynolds and wavenumbers for both of these 
instability modes are the same: Re^ « 642i?a^/^ 8.9 x lO^ffa"^/^ and « 0.477ffa^/^. 
The other pair of modes, which differs from the previous one by an odd distribution of 
vorticity along the magnetic field, is more stable with approximately four times higher 
critical Reynolds number. 



1. Introduction 

Understanding instabilities in magnetohydrodynamic (MHD) flows in ducts is of great 



importance for liquid metal flows in blankets for fusion reactors (|Buhleill2007 ). Blankets 



consist of rectangular ducts in which the liquid metal flows in a high, transverse magnetic 
field of between 5 and 10 T. The aim of these devices is to cool plasma chamber and to 
breed and to remove tritium. This can be assisted by mixing of the flow by turbulence if 
it can be sustained in the presence of a magnetic field. A high magnetic field is known to 
damp turbulence by means of the Lorentz force. At the same time, the magneti c field can 



also a ffect the base v elocity profi le in such a way as to create inflection lines (jKakutani 



1964 ) and even jets (jHuntl 119651 ) thus making the flow more unstable. These two com- 
peting effects balancing each other on a certain length scale result in relatively simple 
asymptotics for the instability threshold. The most dangerous perturbations are usually 
associated with the largest length scale on which the magnetic damping becomes compa- 
rable with the viscous one. This happens in the so-called parallel layers with the relative 



thickness 
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Ha ~^^^ jB raginskii' 'i960*) . Linear stabiHty of these layers for a duct with 



insulating walls ( Shercliff 1953 ) has been considered in a quasi- two-dimensional approxi- 
mation bv lPotheratI 1 2007 ). who found the critical Reynolds number Rcc ~ 4.8x 10^-ffa^^^. 
It is an extremely high, however, typical value for the linear stability of exponential veloc- 
ity profiles, which have th e critical Reynolds nu mber around fifty thousand based on the 
boundary layer thickness (jPrazin fc Reidlll98l[ ). This high threshold is of little practical 
rele vance because the instability for exponential velocity profile is known to be subcrit i- 
cal l|Hockindll975h . This is the case also for the stability of Hartmann layer (|Locklll955l ) , 
which is subcritical too ( Lifshits fc Shtern 1979t Moresco fc Alboussiere 2003 ) with the 
experimentally found Reynolds number for the onset of turbulence in stra ight and an- 
nular ducts of rectangular cross-se ctions being respectively around 225 Ha ( Murgatrovdl 
19531 : iBrouillette fc Lvkoudisl[l967l ) and 380. ffa iMoresco fc Alboussierell2004l ) . Marginal 



turbulent flow states have been observed by IShatrov fc GerbethI tom significantly be- 
low the linear stability threshold in numerical simulations of insulating duct flow subject 
to a transverse magnetic field of moderate strength. 

The stability of MHD flows strongly varies with the electrical conductivity of the duct 
walls. For example, Hunt's flow, which develops in a rectangular duct when the walls per- 
pendicular to the magnetic field are perfectly conducting while the parallel ones are insu- 
lating, has a relatively low linear stability threshold Rc n 91Ha^^^ and Rcr » 112 based 
on th e maximum and average velocities, respectively (jPriede. Aleksandrova fc Molokov 
2010h . The low stability of Hunt's flow is due to two strong jets, which develop in a 
sufficiently strong magnetic fle ld along th e insulating walls and attain a velocity Ha 
relative to that of the core flow (jHuntlll965r ). Although the relative velocity of jets reduces 
as ^ Ha^^^ /c with the increase of the wall conductance ratio c > Ha^^^"^ ( WalkeijllQSl ). 
weak jets with the relative vel ocity 0(1) stifl persist at the parallel waU s also in the limit 
of perfectly conducting duct ( Uflvand 1961 : Chang fc Lundgren 196l[ ). The presence of 
jets with inherent inflection points suggests that this flow may also be highly unstable 
similar to Hunt's flow. It is the aim of the present study to investigate linear stability 
of this flow, which is the last basic MHD duct flow conflguration whose linear stability 
may be not only of theoretical but also of experimental relevance. We first investigate 
the case of square duct and find the high-field asymptotics of the instability threshold 
which are then generalized to arbitrary aspect ratios. 

The paper is organised as follows. The problem is formulated in SJ5] below. In [21 we 
present and discuss numerical results for a square duct in a vertical magnetic field as well 
as for ducts with various aspect ratios in both vertical and horizontal magnetic fields. 
The paper is concluded with a summary and comparison with experimental results in fjH 
Our non-standard vector stream function/vorticity formulation is described in Appendix 

m 



2. Formulation of the problem 

Consider a flow of an incompressible viscous electrically conducting liquid with density 
p, kinematic viscosity u and electrical conductivity a driven by a constant gradient of 
pressure p applied along a straight duct of rectangular cross-section with half-width d 
and half- height h subject to a transverse homogeneous magnetic field B. The walls of 
the duct are assumed to be perfectly electrically conducting and the field may be applied 
across either the width or the height of the duct. 

The velocity distribution of the flow is governed by the Navier-Stokes equation 

dtv + [v ■ V)v = -p-^Vp + vV'^v + p-^f, (2.1) 
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Figure 1. The base flow profile in a rectangular duct with perfectly conducting walls subject 
to a strong vertical magnetic field for Ha — 100. 

where f = j x B is the electromagnetic body force involving the induced electric current 
j , which is governed by the Ohm's law for a moving medium 

j=<j{E + vxB). (2.2) 

The flow is assumed to be sufBciently slow so that the induced magnetic field is neg- 
ligible relative to the imposed one, which supposes the magnetic Reynolds number 
Rm — fiQavod 1, where fiQ is the permeability of vacuum and vq is the character- 
istic velocity of the flow. In addition, we assume that the characteristic time of velocity 
variation is much longer than the magnetic diffusion time = ^oad'^, which allows us 
to use the quasi- stationary app roximation leading to E = — V(/), where </> is the electro- 
static potential ( Robert j 1967t ). The velocity and current satisfy the mass and charge 



conservation V • = V • j = 0. Applying the latter to the Ohm's law (|2.2p yields 

V2(/) = B-u), (2.3) 

where = V x d is vorticity. At the duct walls S, the normal (n) and tangential (t) 
velocity components satisfy the impermeability and no-slip boundary conditions Vn\^ — 
and Vr\g = 0. As the walls are perfectly conducting, the tangential electric current 
vanishes and Ohm's law p.2|) yields (plg = const. 

We employ the Cartesian coordinates with the origin set at the centre of the duct, a;, 
y and z axes directed along its width, height and length, respectively, as shown in figure 
[1] and the velocity defined as v = (u, v, w). The problem admits a purely rectilinear base 
flow with a single velocity component along the duct v — (0, 0, u)(a;, y)) which is shown 
in figure [TJa) for a strong vertical magnetic field. 

In the following, all variables are non-dimensionalised by using the maximum velocity 
Wo and the half-width of the duct d as the velocity and length scales, while the time, 
pressure, magnetic field and electrostatic potential are scaled by d'^/i/, pw^, B = \B\ and 
WodB, respectively. The dimensionless paramerters defining the problem are the Reynolds 
number Re = wod/v, the Hartmann number Ha = dB^ a /{pv) and the aspect ratio 
A — h/d. Note that we use the maximum rather than average velocity as the characteristic 
scale because the stability of this flow, as shown in the following, is determined by the 
former. 

Linear stability of this flow is analysed using the same method as in our previous study 
i Priede. Aleksandrova fc Molokov 20101 ). Since the method is based on a non-standard 



vector stream function formulation, it is briefly outlined in the Appendix. 
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Figure 2. Isolines of the base flow {y > 0) and electric current lines {y < 0) for Ha — 10 {x < 0) 
and Ha = 100 {x > 0) shown in the respective quadrants of duct cross-section (a) and the base 
flow velocity profiles at y = for Ha = 0, 10, 30, 100 (b) in a vertical magnetic field. 



The problem was solved by a spectral collocation method on a Chebyshev-Lobatto grid 
with even number of points defined by 2Nx + 2 and 2Ny + 2 for the x- and y-directions, 
where Nx^y — 35 • • • 60 were used for various combinations of the control parameters 
to achieve accuracy of at least three significant figures. Owing to the double reflection 
symmetry of the base flow with respect to a; = and y — planes, small-amplitude per- 
turbations with different parities in x and y decouple from each other. This results in four 
mutually independent modes, which we classify as (o, o), (o, e), (e, o), and (e, e) according 
to whether the x and y symmetry of i/Jx is odd or even, respecti vely. Our classification 
of mo des co rresponds to the symmetrie s I, II, III, and IV used bv lTatsumi fc Yoshimura 
(1993) and Uhlmann fc Nagata ( 20061) . Thus, four independent problems of different 
symmetries are obtained in one quadrant of the duct cross-section with x Ny internal 
collocation points. The size of matrix for each eigenvalue problem is reduced by a factor 
of 16 in comparison to the original problem. Further detail s and the validation of the 
nume rical method can be found in our previous paper (|Priede. Aleksandrova fc Molokov 
2010h . 



3. Results 

Flow in the presence of a vertical magnetic field induces a transverse current in the 
bulk of the duct which, as seen in the lower part of figure [2l^a), almost directly connects 
to the perfectly conducting walls parallel to the magnetic field. However, a small part 
of the current diverts in the corner regions to connect through the Hartmann walls 
perpendicular to the magnetic field. This makes the density of the transverse current 
and so the resulting electromagnetic force, which opposes the constant driving pressure 
gradient, slightly lower at the parallel walls than in the core region. As a result, weak 
jets form along the parallel walls, where the flow becomes slightly faster than in the core 
of the duct. As seen in figure [JJa) , the formation of jets starts with a velocity minimum 
appearing in the centre of the duct at Ha « 10. With the increase of the magnetic field, 
the velocity in the core becomes almost uniform, while the jets become confined in thin 
layers that develop along the walls parallel to the magnetic field and have a thic kness 
decrea sing as ~ Ha~^^^. In a strong vertical magnetic field, asymptotic solution bv iHund 
yields the velocity maxima located in the mid-plane of the duct at the distance 

S^im5{A/Hay^^ (3.1) 
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from the parallel wall, while the ratio of this velocity to that of the core flow approaches a 
constant 0.809. The latter is seen in figurejSJb) to agree well with our numerical solution. 
An interesting feature of these jets is that the extra flow rate associated with the velocity 
over-shoot above the core velocity balances the flow rate deficit at the wall, where the 
velocity falls below that of the core (|Williamaif 962 ). Thus, the relative contribution of 
these jets to the flow rate is not ~ Ha~^^'^ , as one would expect from simple scaling 
arguments, but a higher-order small quantity Ha~^^^ ^ which is less than the flow rate 
correction due to the Hartmann layers ^ Ha^^ . It is confirmed also by our numerical 
solution, which yields the best fit of the fiow rate for one quarter of a square duct 

Q w 0.809 - 0.4377a-\ (3.2) 

where the leading-order contribution due to the core flow matches the asymptotic solu- 
tion. Note that although t he correction is ~ Ha~^, its coefficient is not equal to that in 
the asymptotic solution bv lHuntI (|l965l ). The difference is because the asymptotic solu- 



tion is obtained for a fixed pressure gradient, whereas numerical solution is for a fixed 
maximum velocity, which has a 0{Ha~^) higher-order correction. Thus, the maximum 
velocity taken as the reference one in this study, results in the same order correction in 
the core velocity and, thus, also in the flow rate. Our choice of the maximum velocity as 
the reference one is motivated by the following results, which show that the instability in 
this fiow is associated with the jets at the parallel walls, which seem inherently unstable 
due to the inflection points in their velocity profiles. Our results can be rescaled to the 
average velocity using relation (|3.2p which becomes sufficiently accurate for Ha > 100. 

We s tart with a square duct, in w hich the fiow without the magnetic field is linearly 
stable ( Tatsumi fc Yoshimural[l990l ). The magnetic field renders the base fiow linearly 



unstable at Ha > 9.6 with respect to a perturbation of symmetry type II. This perturba- 
tion is characterised by the vorticity component along the magnetic field being an odd 
function in the field direction and an even function spanwise. As shown in the following, 
the anti-symmetric distribution along the field results in a strong damping when the field 
strength is increased. The marginal Reynolds number at which the maximum growth rate 
for this mode turns zero (5ft[A] = 0) is plotted in figure [3l^a) against the wavenumber for 
various Hartmann numbers. Besides the marginal Reynolds number, neutrally stable per- 
turbations are characterised by their oscillation frequency uj — 3 [A] and the associated 
phase velocity —uj/k. It is useful to consider the latter relative to the characteristic base 
flow velocity given by Re. This quantity defined as —uj/{Rek) is subsequently referred to 
as the relative phase velocity and shown in figure |3l^b) for mode II. Instability appears 
above the critical Reynolds number Rcc, which is defined by the global minimum on the 
neutral stability curve for the respective Hartmann number. With the increase of Ha, 
the critical Reynolds number for mode II in figure [3ja) first quickly drops to a minimum 
Rcc « 1.1 X 10^ at Ha « 14 and then starts to increase. In some ranges of the Hartmann 
number another minimum appears on the neutral stability curve, which may cause the 
critical mode to jump from the first to the second minimum as the latter becomes the 
global one. This switchover between global minima shows up as a jump in both the crit- 
ical wavenumber and frequency, and as a break point in the dependence of the critical 
Reynolds number on the Hartmann number. Such a jump is noticeable in figure [5] at 
Ha w 48, where the instability switches from mode Ha to lib. 

However, this jump is of secondary importance because a mode of type I is seen in 
figure [SJa) to become more unstable than mode II at Ha « 10.6. Although mode I 
turns linearly unstable at a slightly higher Hartmann number than mode II, its critical 
Reynolds number decreases faster with the increase in the Hartmann number than that 
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Figure 3. The marginal Reynolds number (a) and the relative phase velocity (b) versus the 
wavenumber for neutrally stable modes of type II in a square duct (^4 = 1) subject to a vertical 
magnetic field at various Hartmann numbers. 



for mode II. As seen from the neutral stability curves in figure |3l^a,b), with the increase 
of the Hartmann number, the critical Reynolds number for mode I first quickly drops 
to to a minimum Rcc « 2697 at Ha w 17 and then starts to raise. With the increase 
of Ha mode I is quickly approached from above by a mode of type III, which is seen 
in figure EJa) to appear at Ha 30 and become practically indistinguishable from the 
former at Ha > 80. There is also a mode of type IV appearing at Ha w 30, which in 
turn approaches mode II in a similar way as mode III approaches mode I. Modes III/IV 
differ from modes I/II by the opposite symmetry across the magnetic field. Namely, for 
modes III/IV, the vorticity component along the magnetic field is an odd function in the 
spanwise direction across the magnetic field, whereas it is an even function for modes 

i/n. 

The most important feature of the instability seen in figure [5]^ a, 6) is the critical 
Reynolds number and the wavenumber for each of two merged pairs of modes increas- 
ing in strong magnetic field as ^ Ha^^'^. The relative phase velocity shown in figure 
[5] tends to a constant ^ 0.911 for both pairs of modes. This kind of variation implies 
that in a strong magnetic field the instability is determined by the internal length scale 
6 ~ Ha~^^^ , which is the characteristic thickness of the jets developing along the walls 
parallel to the magnetic field. The base flow velocity correction of order Ha^^, which 
was discussed above, implies an 0{Ha~^^^) correction to the critical Reynolds number. 
The best fit for modes I/III yields 

Rec{,Ha\ A = 1) w M2Ha^/^ + 8.9 x 10^7^a"^/^ (3.3) 
kc{Ha; A = 1) w 0.4777^a^/^ (3.4) 

which are seen in figure [5j a, &) to well approximate numerical results for mode I down to 
Ha « 30. Similarly, for modes II/IV, we find 

Re^{Ha- A^l)^ 2580ifa^/^ + 1.1 x IQ^Ha^^'^, (3.5) 
kc{Ha] A = 1) w 0.4197^a^/^ (3.6) 

where the former is by nearly of a factor of four greater than (|3.3p . 

The instability threshold being nearly the same for the modes of opposite spanwise 
symmetry in a sufficiently strong magnetic field implies that the perturbations developing 
in the jets at the opposite walls are effectively separated by the core region of the flow 
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Figure 4. The marginal Reynolds number (a,c) and the relative phase velocity {b,d) versus 
the wavenumber for neutrally stable modes of type I {a,b) and type III (c, rf). 



and, thus, do not affect each other. This is confirmed by the patterns of the critical 
perturbations, which are plotted over the duct cross-section in figure [S] for a moderate 
{Ha = 15) and a relativel y strong (Ha = 100) magnetic field. As shown in our previous 
paper ( Priede. Ale ksandro va fc Molokovll2010( ). the flow perturbation can be represented 



by the complex amplitudes of the streamwise (z) component of velocity (w) and that of 
stream function (ipz), whose isolines are plotted in the left {x < 0) and the right (a; > 0) 
sides of the cross-section, respectively. The real and imaginary parts of perturbations 
plotted at the top and bottom halves of the cross-section show the instant patterns shifted 
in time or in the stream- wise direction by a quarter of period or wavelength, respectively. 
Although the real and imaginary parts of the complex amplitude distributions completely 
determine the evolution of perturbation over the harmonic oscillation cycle, these two 
quantities are not uniquely defined. The main ambiguity is due to the free choice of the 
initial time instant and the initial stream-wise coordinate. This ambiguity can be partly 
eliminated by choosing the phase of the complex velocity perturbation amplitude so that 



where the integrals are taken over the duct cross-section S. This condition definines the 
phase up ±7r/2, which means swaping the real and imaginary parts. The perturbation 
amplitude remains defined up to a constant factor which is not important as the linear 
stability theory predicts only the pattern but not the amplitude of the critical perturba- 
tions. 
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Figure 5. The critical Reynolds number (a), wavenumber (b) and relative phase velocity (c) 

against the Hartmann number. 



Positive and negative values of w are respectively associated with converging and di- 
verging potential flow component in the cross-section plane. The isolines of i/j^ correspond 
to the streamlines of the solenoidal flow component in that plane. The critical perturba- 
tions for modes I and II, which are shown in figures a, &) and {d,e), respectively, differ 
by their vertical symmetry. Namely, the perturbation of w and i/'^ are respectively even 
and odd functions of y for mode I, whereas they are odd and even functions for mode II. 
Thus, the vortices for mode I rotate in opposite senses in the upper and lower parts of 
the cross-section, whereas for mode II there is one symmetric vortex spanning the whole 
height of the duct. For both of these modes, the pairs of vortices across the vertical mid- 
plane rotate in the same sense and, thus, represent two parts of a bigger vortex spanning 
over the whole width of the duct. At Ha = 15, slightly above the Hartmann number 
at which the flow turns linearly unstable, the critical perturbations are seen in figure 
[n{a,(f) to be localised close to the duct centre, where the two velocity maxima discussed 
above first appear. In this case, the co-rotating vortices on the opposite sides of the duct, 
whose symmetric half is shown at a; > 0, are clearly connected by the flow through the 
vertical mid-plane. However, this is no longer the case for a sufficiently strong magnetic 
field. As seen in figure in[Sl^&,e), at Ha = 100 the critical perturbations for modes I and 
II are localised at the side walls. Moreover, these perturbations have nearly the same 
pattern as those for modes III and IV, which are shown in figure [6jc,/) for the same 
Ha. It means that in a strong enough magnetic field the perturbations at the opposite 
walls are effectively separated by a stagnant core flow and, thus, do not affect each other. 
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Figure 6. Amplitude distributions of the real {y > 0) and imaginary [y < 0) parts of w (x < 0) 
and (x > 0) of the critical perturbations over one quadrant of duct cross-section for instability 
modes I {a,b), Ila (d,e), III(c) and IV(/) at Ha= 15 {a,d) and Ha = 100 {b,c,e,f). 



This explains the merging of the instability thresholds for the modes with the opposite 
spanwise symmetries seen in figure [Sj 

In weaker magnetic fields or in narrower ducts across the magnetic field, which will be 
considered later, perturbations at the opposite walls can either enhance or suppress each 
other depending on their spanwise symmetry. The first is the case for the perturbations 
of type I/II, which, as discussed above, have co- rotating vortices at the opposite walls 
connected by a fiow across the vertical mid-plane. These perturbations are more unstable 
than those of type II/IV with counter-rotating vortices at the opposite walls, which tend 
to suppress each other, especially in moderate magnetic fields or in sufficiently narrow 
ducts. 

Besides the spatial amplitude distributions perturbations can be characterised quan- 
titatively by their kinetic energy distribution, which can be represented in two different 
forms using either the velocity or vorticity /stream function components 



Js Js 

where E is the kinetic energy of pertu rbation averaged over the wavelength and the 
asterisk denotes the complex conjugate ( Priede. Aleksandrova fc Molokov 20101 ). At the 
moderate Hartmann number Ha = 15 considered above, most of the kinetic energy, i.e. 
60% and 66% for mode I and 88% and 53% for mode II, is carried by the z-component 
of velocity and by the y-component of vorticity, respectively. For mode I. next most 
energetic is the x-component of both velocity and vorticity, which contain respectively 
about 27% and 23% of the energy. For mode II, in this position are the y velocity and the 
X vorticity components, which contain respectively about 10% and 37% of the energy. 

The energy distribution for the most unstable modes I and II becomes simpler in a 
strong magnetic field. For example, at Ha = 100, 86% of the energy for these modes is 
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Figure 7. Isosurfaces ofipy (a) and longitudinal velocity w (b) perturbations over wavelength in 
one quadrant of the duct cross-section for instability mode I at Ha — 100 in a vertical magnetic 
field. 

concentrated in the y-component of vorticity, while the rest is distributed nearly equally 
between the two other vorticity components. This component of vorticity is associated 
with the circulation in the (x, z)-planes transverse to the magnetic field. The streamlines 
of this circulation are represented by the isolines of t/jj,, whose spatial distribution is 
shown in figure [71[ a). Note that tlie distribution ofipy is very close to that of the electric 
potential </> because the equations (jA 6[) and (|A 5[) governing both quantities are the 
same. Moreover, ^py and </> satisfy the same boundary condition at the wall parallel to 
the magnetic field. Thus, both quantities differ only in the vicinity of the wall normal to 
the magnetic field, where they have different boundary conditions. 

The transverse character of circulation for modes I and III is confirmed also by the 
kinetic energy distribution over the velocity components. Only about 9% of the energy is 
carried by the velocity component along the magnetic field, while 67% by the streamwise 
(z) velocity perturbation, whose spatial distribution is shown in [TJ^b) . The relatively low 
contribution of the spanwise (x) velocity component, which carries the remaining 24% of 
the energy, is due to the relatively long wavelength of perturbation Ac = 27r/fcc, which 
according to (|3.4p is by approximately a factor of 7 larger than the thickness of the jet 

(EH). 

In a strong magnetic field, the energy distribution in modes II/IV is essentially differ- 
ent from that in modes I/II considered above. Namely, the latter two have only 5 — 6% 
of their energy in the spanwise (x) velocity component, which implies a circulation con- 
strained mainly to the (y, z)-planes parallel to the side walls. Similar to the previous two 
modes, 66% of the energy is carried by the streamwise velocity component. Although 
circulation mostly occurs in the (y, 2;)-planes, only 19% of the energy is contained the 
transverse (x) vorticity /stream function component, while 56% are still contained in the 
vorticity /stream function component along the magnetic field. This scatter of energy 
between the vorticity components is due to the confinement of circulation in narrow 
layers parallel to the side walls. The confinement causes a strong variation of the veloc- 
ity perturbation in the spanwise (x) direction and, thus, produces a significant vorticity 
components tangential to the plane of circulation. 

Finally, we consider the effect of the duct aspect ratio on the instability threshold of 
the two most unstable modes. In order to generalise the above results for square duct to 
arbitrary aspect ratios it is instructive to start with a horizontal magnetic field, which is 



Linear stability of MHD flow in a perfectly conducting rectangular duct 1 1 

directed along the fixed dimension of the duct, i.e. the width in our case. In sufficiently 
strong horizontal magnetic field, the variation of the instability threshold with the aspect 
ratio turns out to be particularly simple and directly related to that of the square duct 
considered above. We shall use this fact and the similarity of horizontal and vertical 
magnetic field orientations to obtain more general relations for the latter orientation. 

It is important to notice that the thickness of parallel layers 5 defined by p.ip . which 
is the characteristic length scale of the instability, varies not only with the Hartmann 
number but also with the aspect ratio A, which defines the size of the duct along the 
magnetic field. When the magnetic field is directed horizontally along the fixed dimension 
of the duct, 5 becomes independent of A and varies only with Ha as in the case of square 
duct. This simplification is our main motivation for considering first horizontal magnetic 
field. 

On changing the magnetic field from vertical to horizontal, modes I/III swap with 
II/IV, which, thus, become the most unstable ones. For sufficiently large aspect ratios, 
the critical Reynolds number and the wavenumber are seen in figure [51^ a, 6) be the same 
for both modes, which agree with the strong field asymptotics p.3l3.4p . As discussed 
above, this implies that the instabilities developing in the jets at the parallel walls are 
effectively separated by a stagnant core of the flow and, thus, do not affect each other. It 
changes at small aspect ratios, when the walls parallel to the magnetic field are sufficiently 
close to each other. Then the vortices at the opposite walls start to interact, which causes 
the thresholds for both modes to diverge. For mode IV, the vortices at the opposite walls 
counter-rotate and, thus, tend to suppress each other, which results in the stabilisation 
of the flow. It is the opposite for mode II, whose instability threshold flrst drops as the 
co-rotating vortices at the opposite walls start to enhance each other. However, with 
further reduction of the aspect ratio the critical Reynolds number attains a minimum 
and then starts to increase following that for mode IV. The raise of Rcc for mode II is 
associated with the increase of the critical wavenumber. This corresponds to the reduction 
of wavelength which is required for the vortices to fit between closely spaced parallel walls. 

Now we turn to vertical magnetic field, which is oriented along the variable height of 
the duct. This slightly more complicated case can be reduced to the previous one by 
taking the height of the duct as the length scale and rotating the duct by 90 degrees. 
This is equivalent to the substitutions A' = 1/A, Re — ARe, Ha' — AHa, and k' — Ak, 
where the parameters with primes correspond to the case of horizontal magnetic field 
considered above. Then we can use the result for horizontal magnetic field found above 
according to which the critical parameters in strong magnetic field approach those for 
a square duct given by p.3l3.4p . Substituting the primed parameters into p.3l3.4p we 
obtain 

RcciHa; A) = A-^RedAHa; 1), (3.7) 
kciHa; A) = A'^kdAHa; 1), (3.8) 

which are seen in figure [8l^c,d) to fit the numerical results well in the intermediate range 
of aspect ratios, where the thresholds for both most unstable modes I/III merge. In this 
range both the critical Reynolds number and wavenumber for a fixed Ha vary asymptoti- 
cally with the aspect ratio as ~ A~^/'^. This variation, which is due to the increase of the 
jet width p.ip as ^ A^/^, breaks down at both small and large A. In the former limit, 
asymptotic relations (|3.7l3.8p turn inapplicable because the effective Hartmann number 
AHa based on the height becomes too small for (|3.3|3.4|) to be valid. In the latter limit, 
the jets become so wide that the vortices at the opposite walls start to interact resulting 
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(a) (b) 




Aspect ratio, A Aspect ratio, A 

Figure 8. The critical Reynolds number (a,c) and the wavenumber (b,d) versus the aspect ratio 
for modes II/IV in a horizontal magnetic field (a,b) and for modes I/III in a vertical magnetic 
field (c,d) at various Hartmann numbers. 

in the divergence of the instability thresholds and the eventual stabilisation described 
above for the case of horizontal magnetic field. 

4. Discussion and conclusions 

We have presented numerical results concerning linear stability of a liquid metal flow 
in a rectangular duct with perfectly electrically conducting walls subject to a uniform 
transverse magnetic field. It was found that a linearly stable flow in a square duct turns 
unstable as a relatively weak magnetic field with the Hartmann number Ha « 9.6 is 
applied. The instability is due to two weak jets, which first appear at the centre of the 
duct and then move to the walls parallel to the magnetic field as the field strength is 
increased. The instability follows the jets and in a sufficiently strong magnetic field be- 
comes confined in the layers of characteristic thickness S ^ Ha~^^^ located at the parallel 
walls. The thickness S determines the characteristic length scale of the instability, which 
results in both the critical Reynolds and wave numbers scaling as ~ 5~^. Owing to the 
double reflection symmetry of the problem, perturbations with four different symmetry 
combinations along and across the magnetic ffeld decouple from each other and, thus, 
are considered separately. The most unstable is a pair of perturbations with an even dis- 
tribution of the vorticity along the magnetic ffeld. These two modes represent strongly 
non-uniform vortices aligned with the magnetic field, which rotate either in the same or 
opposite directions across the magnetic field. The former enhance while the latter weaken 
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one another provided that the magnetic field is not too strong or the walls parallel to the 
field are not too far apart. In a strong magnetic field, when the vortices at the opposite 
walls are well separated from the core flow, the critical Reynolds number and wavenum- 
bers for both of these instability modes are the same: 

1/2 

kr ~ Q.mHa''\ The other pair of critical perturbations, which differs from the previ- 
ous one by an odd distribution of vorticity along the magnetic field, is more stable with 
approximately four times higher critical Reynolds number. 

The basic instability characteristics described above resemble those for the Hunt's 
flow, which has a simil arly increasing, however a significantly lower, critical Reynolds 
number Rcc ^ 91Ha^^^ i Priede. Aleksandrova &: Molokov 2010t ). The difference becomes 
substantial when the average rather than the maximum velocity is considered. Namely, 
the critical Reynolds number based on the average velocity for this flow increases in a 
strong magnetic field in the same way as Rcc ~ 519Ha^^^ , whereas it tends to a constant 
Rec w 112 for Hunt's flow. This asymptotically constant Rcc is due to a principally differ- 
ent velocity distribution in Hunt's flow with the jets o f thickness ^ Ha~ ^^^ dominating 
the flow rate. Constant Rcc ~ 313 has been found bv iTing et al. (1991) for the linear 
stability of the flow in a square duct with thin but relatively well-conducting walls in 
a strong transverse magnetic field. This flow represents an intermediate case in terms 
of the conductivity of parallel wall between the perfectly conducting one considered in 
this study and the insulati ng one for Hunt's flow. Although the leading-order velocity 
distribution considered by iTing et al\ (|l99ll ) is nearly the same as that of Hunt's jet, 
there is one principal difference between two. Namely, in a duct with thin parallel walls 
of a moderate conductance ratio c — Uwdw/crd, where and dw are the electrical con- 
ductivity and the thickness of wall, satisfying Hd~^^^ <^ c <^ Ha}^^^ ^ jets carry a volume 
flux of the same order of magnitude as that of the core flow, whereas the contribution 
of the latter is negligible in Hunt's flow. In a square duct with thin walls, in which the 
core flow in strong magnetic field carries 3/4 of the total volume flux (Walker 1981), the 
jet velocity is by a factor of about four lower than that for Hunt's flow at the same Re. 
It means that jets are more unstable at wall of finite conductivity than in both limiting 
cases of insulating and perfectly conducting parallel walls. 

The experiment best matching the model considered in this study has been carried 
out bv lBranover fc Gelfgat ( 1968 ). who measured a flow of mercury at i?e « 4 x 10''' and 
Ha = VIA in a rectangular copper duct with a relatively high wall conductance ratio, 
which was 20 and 10 for the Hartmann and the parallel walls, respectively. The magnetic 
field was applied transversely to the longest side of the duct with the aspect ratio of 1.5. 
First, the authors found a maximum velocity in jets exceeding the theoretical prediction 
for perfectly conducting duct by approximately 50%. At this Hartmann number, we find 
that 19% of the excess jet velocity may be due to the finite wall conductivity, when the 
latter is included in the numerical solution. Note that the effect of imperfectly conducting 
walls is not entirely det ermined by t he relative conductance of the par allel layers cHo}^'^ 
as it may appear from Hunt ( 1965 ) (see also Miiller fc Biihler 200ll p. 146). Namely, 
cHa^^^ ^ 1 means only that the parallel walls cannot be treated as insulating (c ^ 0). 
But it does not necess ary mean tha t the walls may be assumed perfectly conducting 
(c = cxd). As shown bv IWalkei (1981), the latter approximation requires a much higher 
wall conductance ratio cHa ' ^1. Thus, the effect of imperfectly conducting walls 
increases rather than decreases with the field strength, which results in t he jet velocit y 
relative to that of the core increasing with the field strength as Ha^^^ /c ( Walkeij [ l981 ). 
However, the most significant deviation f rom the laminar flow solut ion by Hunt! ( "9651 ) 
was the jet thickness, which was found by lBranover fc Gelfgat ( 1968h to be several times 
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greater than expected. Such a broadening of jets is Ukely to be due to the turbulence 
which is expected in the experiment at the Reynolds number significantly above the 
linear stability threshold Rcc ~ 6850 predicted by our analysis for this setup. 

In conclusion, note that when the distance of the velocity maximum from the parallel 
wall (|3.ip is taken as the length scale, p.3p becomes i?ef w 1230. This small critical 
Reynolds number may be due to the inviscid nature of the instability caused by the 
inflection point of the base velocity profile. It also implies that the instability may be su- 
percritical rather than subcritical as for typical shear flows. The supercritical character of 
instability may expl ain why it appears significant l y abov e the linear stability threshold as 
recently reported bv lKinet. Knaepen fc Molokovl ( 20091 ). They observed small- amplitude 



vortices in the jets at the parallel walls for 2500 < Re < 3700 at Ha — 200 in the numer- 
ical simulation of a flow in a rectangular duct with thin walls. These vortices were found 
to be subcritically unstable for 3500 < Re < 3700, whi ch may co rrespond t o the large 
ampU tude instabilities observed in the experiments (Ree d fc Picolog lou 19 89[ , lBurr etal 



2000). If the instability in perfectly conducting duct like that in ducts with thin walls of 



finite conductivity is supercritical, it is expected to appear above the absolute instability 
threshold which cannot be lower than that of the co nvective instability pred icted by the 
classical linear stability analysis used in this study ( Priede fc Gerbeth IQQTI) . The abso- 



lute instability threshold for these jet-type flows, at which small-amplitude self-sustained 
vortices are expected to appear, is not yet known. 
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the opportunity to use its high performance computer cluster. 



Appendix A. Vector stream function/ vorticity formulation 

We use the vector stream function xp, which is introduced to satisfy the incompressiblity 
constraint V • v = for the flow perturbation by seeking the velocity distribution in the 
form V = V X '0. Since the velocity is invariant upon adding the gradient of arbitrary 
function to xp, we can impose an additional constraint 

V-ip = 0, (Al) 



which is analogous to the Coulomb gauge for the magnetic vector potential A i JacksonI 



IQQSl ). Similarly to the incompressiblity constraint for v, this gauge leaves only two 



independent components of xp. 

The pressure gradient is eliminated by applying curl to (|2.ip . This yields two dimen- 
sionless equations for xp and lo 

dtu: ^V^u} ~ Reg + Ha^h, (A2) 
O^V^xp + u;, (A3) 

where g — 'V x {v ■ V)v, and /i = V x / are the curls of the dimensionless convective 
inertial and electromagnetic forces, respectively. 

The boundary conditions for xp and u are obtained as follows. The impermeability 
condition applied integrally as j^v ■ ds = §^xp ■ dl — Q to an arbitrary area of wall s 
encircled by a contour I yields iPt\s =0. This boundary condition substituted into (jA 1|) 
results in dn''pn\g = 0. In addition, the no-slip condition applied integrally f^v ■ dl = 
JgU: ■ ds yields a;„|^ =0. 

Linear stability of the base flow {xp, u>, (/>}(a;, y) is analysed with respect to infinitesimal 
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disturbances in the form of harmonic waves travehing along the axis of the duct 

{tA, a;, 0}(r, t) = {i,, u,, 0}(x, y) + iJ, $}{x, y)e'^*+«'=^ 

where A: is a wavenumber and 7 is, in general, a complex growth rate. This expression 
substituted into (jA2|A3P results in 

-ids^VlCj- Reg + Ha^h, (A 4) 

= V|i^^ + w, (A 5) 

O = Vl0-u;,„ (A 6) 

where V/c = V^+fi/cez; 11 and ± respectively denote the components along and transverse 
to the magnetic field in the (a;, y)-plane. Because of the solenoidality of cIi, we need only 
the X- and y-components of (jA4|) . which contain h±^ ~ ~dxy(f> — d„w, h„ — —d'^(j> and 

§3: — k^vw + dyy{vw) + dxy{uw) + &2kdy{ww) , (A 7) 

cjy = —k^uw — dxx{uw) — dxy{vw) — Wlkdxiww), (A 8) 

where 



M = 6fc ^{dyylpy - k'^ljjy + dxylPx), 

V = -?>k^^{dxxi^x - k^-ipx + dxyipy), 



W = dxtjjy - dylpx- 

The relevant boundary conditions are 

(j)^^y^ dxi^x dxtpy ~ dyijjx = cDj, = at a; = ±1, (A 9) 

<^ = V'x = dy->jjy = dxipy - dy'4)x = cjy = at J/ = ±A. (A 10) 
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